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We discuss various mechanisms for the creation of an asymmetric charge fluctuation with respect 
to the reaction plane among hadrons emitted in relativistic heavy-ion collisions. We show that such 
mechanisms exist in both, the hadronic gas and the partonic phases of QCD. The mechanisms con- 
sidered here all require the presence of a strong magnetic field (the "chiral magnetic effect"), but 
they do not involve parity or charge-parity violations. We analyze how a transient local electric cur- 
rent fluctuation generated by the chiral magnetic effect can dynamically evolve into an asymmetric 
charge distribution among the final-state hadrons in momentum space. We estimate the magni- 
tude of the event-by-event fluctuations of the final-state charge asymmetry due to the partonic and 
hadronic mechanisms. 



I. INTRODUCTION 

Collisions of two heavy nuclei at high energy serve 
as a means for creating and exploring strongly inter- 
acting matter at the highest possible energy densities 
where matter is expected to assume the state of a quark- 
gluon plasma [l| . The properties of matter governed by 
the laws of quantum chromodynamics (QCD) have been 
studied in this way for a decade at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory 
The measurements performed in Au+Au, Cu+Cu, 
d-|-Au and p+p collisions at centcr-of-mass energies up 
to 200 GcV per nuclcon pair have revealed several un- 
usual properties of such supcrdense, strongly interacting 
matter [6|, most notably its very low kinematic shear 
viscosity and its high opacity with respect to energetic 
particles carrying a free color charge. 

The strong interactions are known to respect space and 
time reflection symmetry to a very high degree. This is 
not a direct consequence of the laws of quantum chromo- 
dynamics which, in principle, permit a so-called 0-tcrm 



which violates time reversal symmetry. (Here stands 
for the gluon field strength tensor and F^a"' for its dual.) 
Rather, the symmetry conserving nature of QCD has 
been established by precise experiments that set limits 
on the intrinsic electric dipole moment of the neutron. 
The present experimental limit implies that the co- 
efficient of the possible CP violating term in the QCD 
Lagrangian 6 < 0.7 x 10~^^ Q. The reason for its sup- 
pression or even its complete absence is not known; an 
often considered mechanism is the postulated existence 
of a new, spontaneously broken symmetry, called Peccei- 
Quinn symmetry which would give rise to a new light, 
neutral particle, the axion [loj . 

Even if CP-symmetry is not violated in the normal 
QCD vacuum, it is conceivable that it is violated in an 
excited state or "false" QCD vacuum. Theorists have 



speculated that such a state might be created when QCD 
matter is heated above the critical temperature for con- 
finement and chiral symmetry breaking and accidentally 
cools down into a configuration that is characterized by a 
nonzero value of 9 [111, Illl , and experimental observables 
capable of detecting the formation of symmetry violating 
domains ofQCD matter in heavy ion collisions have been 
proposed |13l4l5l |. 

Colliding nuclei create not only the highest tempera- 
tures, but also the strongest magnetic fields attainable 
in the laboratory [l6j . For the purpose of studying the 
properties of QCD matter, the coherent magnetic field 
generated by the moving nuclei can be considered as "ex- 
ternal" : the value of the magnetic field at a given point 
is determined by the global charge distribution of the 
colliding nuclei and is thus, in good approximation, inde- 
pendent of the local strong interaction dynamics. How- 
ever, the presence of the magnetic field can affect the 
interactions among the quarks and anti-quarks, which si- 
multaneously carry electric and color charge. 

Because magnetic fields are odd under time reversal (or 
cquivalently, under a combined charge conjugation and 
parity (CP) transformation), the time reversal symme- 
try of a quantum system is broken in the presence of an 
external magnetic field. A magnetic field B can also com- 
bine with an electric field E to form the Lorentz invariant 
E • B which changes sign under a parity transformation. 
An especially interesting aspect of magnetic fields is thus 
that they can be used to probe the response of strongly 
interacting matter in the pseudoscalar, CP-odd sector of 
QCD. In the simultaneous presence of parallel electric 
and magnetic fields, QCD matter is thus dynamically 
forced to explore states with unnatural CP symmetry, 
eliminating the need to invoke a "spontaneous" symme- 
try violating transition. 

As first pointed out by Kharzeev, McLerran, and War- 
ringa [2l|, the coherent magnetic field generated by two 
heavy nuclei colliding off-centrally at high energy can 
convert topological charge fluctuations in the quantum 
chromodynamic (QCD) vacuum into global fluctuations 
around the electric charge symmetry with respect to the 
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reaction plane. They called this mechanism the chi- 
ral magnetic effect and interpreted it "local" vio- 

lation of P and CP symmetry. The chiral magnetic effect 
was further analyzed by Fukushima, Kharzeev, and War- 
ringa [23 - [23 |. who argued that it is specific to the high 
temperature phase of QCD matter, commonly called the 
quark-gluon plasma, in which quarks arc liberated from 
confinement into hadrons, because only then could the lo- 
cal electric current fluctuations be converted into charge 
fluctuations on the nuclear scale. The same conclusion 
was reached by Fu, Liu, and Wu , who calculated the 
electric charge separation in a magnetic field in the frame- 
work of the Polyakov loop-Nambu-Jona-Lasinio (PNJL) 
model. The temperature dependence of the chiral mag- 
netic effect and its dependence on the magnetic field in 
the presence of instanton fiuctuations was investigated 
by Nam [2^. Lattice QCD simulations of the chiral mag- 
netic effect were reported by Buividovich et al. [23] (in 
quenched SU(2) gauge theory) and Abramczyk et al. [2^ 
(in (2-|-l)-fiavor QED+QCD). A review of the chiral mag- 
netic effect and its possible manifestation in relativistic 
heavy- ion collisions can be found in [2^ . 

Inspired by these theoretical arguments and by exper- 
imental considerations related to the detection of parity 
violating effects in QCD [sO], the STAR collaboration has 
analyzed the final states of heavy-ion collisions at RHIC 
(Au+Au and Cu+Cu) for the presence of non-vanishing 
fluctuations in P- and CP-odd observables. The most 
promising of these is the asymmetry with respect to the 
reaction plane of the average emission angle between like- 
sign and opposite-sign hadrons defined by the event av- 
erage 

C(±'^' = (cos(0W+0(±)-2*Hp)) (2) 

where (f>a, 4>p denote the azimuthal emission angles of any 
pair of hadrons with respect to the beam axis, and 4'rp 
denotes the azimuthal orientation of the reaction plane. 
A difference 

AQ - + ' - 2C(+-) (3) 

indicates a charge asymmetry of the kind suggested by 
the chiral magnetic effect. The existence of such an asym- 
metry was reported in a recent publication [3l[ . 

The purpose of our manuscript is to address two is- 
sues not clearly discussed in previous publications. First, 
we wish to point out that the mechanism producing an 
asymmetric charge fluctuation of the kind induced by the 
chiral magnetic effect is not constrained to an environ- 
ment in which quarks are dcconfined. We show that an 
analogous mechanism, magnetic tt — p conversion, exists 
in a confined hadronic environment. This should be ex- 
pected in view of the general principle of parton-hadron 
duality, which is pervasive in QCD |32l | and asserts that 
every mechanism present at the partonic level has an 
analogue in the world of hadrons. 

A second purpose of our manuscript is to analyze how 
a local electric charge asymmetry can be converted into 



a global asymmetry of charged hadron emission. This is 
a nontrivial problem, because the electric current fluctu- 
ations with respect to the reaction plane constitute an 
effect in position space, whereas the observable ([5]) mea- 
sures an asymmetry in momentum space. We show that 
the presence of transverse collective flow and, more gen- 
erally, the properties of hadronic freeze-out are essential 
to converting an asymmetry in position space into a mo- 
mentum space asymmetry. 

We begin by discussing the partonic and hadronic 
mechanisms for the electric current fluctuations under- 
lying the chiral magnetic effect (CME) and then investi- 
gate the question of how a fluctuation in the electric cur- 
rent with respect to reaction plane gets converted into a 
fluctuation of the emission of charge hadrons. Next, we 
estimate the size of certain hadronic and partonic (glu- 
onic) contributions to a charge asymmetry fluctuation. 
Our work concludes with a summary and outlook. 



II. PSEUDOSCALAR INTERACTIONS OF THE 
ELECTROMAGNETIC FIELD 

A. Role of the Axial Anomaly 

In the normal QCD vacuum with its spontaneously 
broken chiral symmetry, the leading interaction involving 
the invariant E ■ B is (l7l.[l8j: 

C= ^-/-^E-B, (4) 

t—7T^' ,7] ,7] 

where (pi denotes one of the neutral pscudoscalar me- 
son fields, fi is the respective meson decay constant, and 
a stands for the electromagnetic fine-structure constant. 
The interaction ^ mediates the two-photon decays of 
the neutral pscudoscalar mesons. 

In the deconfined, chirally symmetric phase of QCD, 
the leading interaction is of the form 

£' = Kaa,(E'' •B°)(E-B), (5) 

where E" and B" denote the chromoelectric and chro- 
momagnetic fields, respectively, and as = g^/47r is the 
QCD coupling constant. The interaction ^ is not fun- 
damental to QCD; it arises as an effective interaction at 
the one-quark- loop level (lOj . 

Both interactions are closely related to the electromag- 
netic axial anomaly, which relates the divergence of the 
isovector axial current to the pscudoscalar invariant of 
the electromagnetic field 

'^mjF'"' = i{'^uU-f5U + nidd-fsd) - ^^^^^E • B, (6) 

where m, are the current quark masses and N^, = 3 is 
the number of colors. The fiavor-singlet (isoscalar) axial 
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current in QCD has the anomaly 
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Nfa. 



2Nca 



(7) 



where Nf the number of hght flavors, e/ denotes the 
electric charge of each quark flavor, and ag = g'^/in is 
the QCD coupling constant. (For a recent review of the 
physics of the axial anomaly see [l^l-) The pseudoscalar 
invariant of the gluon field is related to the topological 
charge density (or Chern-Simons number density) of the 
gluon field 
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(8) 



B. Anomalous Current in QCD 

The "chiral magnetic effect" occurs when the dynam- 
ics governing the nuclear reaction induces an electromag- 
netic current fluctuation 5] that is parallel to the mag- 
netic field B generated by the colliding nuclei. The elec- 
tromagnetic current operator is given by 



(9) 



where £off[^] is the effective Lagrangian for the electro- 
magnetic field. In order to induce a current parallel to 
the external magnetic field B, £cff must contain a con- 
tribution of the pseudoscalar form 



(10) 



where F'"' = (l/2)e'"''"'3Fa/3 is the dual field tensor and 
P is a pseudoscalar operator constructed from strongly 
interacting fields. Using this form of the effective La- 
grangian, one obtains the following general form for the 
anomalous electromagnetic current: 



,(x) = -(9,7')F^'^ 



(11) 



which was first derived in general form by D'Hoker and 
Goldstone [s^. By virtue of Maxwell's equations and 
the antisymmetry of F'"^, one easily confirms that the 
anomalous current is conserved: 

d^.f.A^) = ~{d^.d.V) F^"- - {d,V) d.F^'' = 0. (12) 




FIG. 1: Feynman diagram describing the anomalous coupling 
between two photons and a neutral pseudoscalar meson via a 
triangular quark loop. 



pseudoscalar mesons (i = 7r°, 77,77'), which give rise to 
interactions of the form ([TU]) via the Adler-Bell-Jackiw 
triangle anomaly (see Fig. [1]): 



(had) 
eff 



E^^^(E-B), 



(13) 



where (j^i denotes the pseudoscalar meson fields and fi 
stands for the meson decay constant. In the chirally 
unbroken, deconfined phase of QCD the relevant pseu- 
doscalar operator is given by the topological charge (or 
Chern-Simons number) density pes defined in ([8]). The 
combined effective QED-QCD action for soft gauge fields 
contains a term of the form 1191 



(QECD) 
off 



Kaas (E'^-B'") (E • B). 



(14) 



The coupling coefficient n oc X]/(^//^)^' where e/ de- 
notes the electric charge of a quark of fiavor /, is given 
by the quark box diagram (see Fig. [2]) with two gluon 
and two photon vertices. 




FIG. 2: Feynman diagram describing the effective pseu- 
doscalar coupling between two gluons and two photons via 
a quark loop. 



As mentioned in the Introduction, the form of possi- 
ble pseudoscalar operators V differs fundamentally be- 
tween the phases of QCD with and without spontaneous 
chiral symmetry breaking. In the normal QCD vac- 
uum with spontaneously broken chiral symmetry, rele- 
vant pseudoscalar operators exist in the form of neutral 



In the limit that the magnetic field, the current, and 
the exchanged gluons are "soft" and the quark mass tti is 
large, the quark box diagram is obtained from the gen- 
eralization of the Heisenberg-Euler effective Lagrangian 
of QED to include both, electromagnetic and Yang-Mills 
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fields. At fourth order in the fields, the eff'ective La- 
grangian in the vacuum reads (l9j : 



(4) 



-1 



Tr 



^QECD - 3g07,2 

where 

f- — 1 F F^^'^ 



g = 



(16) 



are the two Lorcntz invariants of the "mixed" field tensor 




FIG. 3: Feynman diagram describing the effective pseu- 
doscalar coupling between two gluons and two photons via 
a virtual 77- or 77'-meson. 



F = eF 



(17) 



Here F^j^y denotes the electromagnetic field tensor, F^^ 
denotes the field tensor of the color field, and the Gell- 
Mann matrices t"" are generators of color-SU(3) in the 
fundamental representation. The trace in (|15p runs 
over colors and quark flavors; the relevant color trace 
is TT{tH^) = (5°V2- The contribution of interest to us is 



Trc(^') = ^e'g\F;,F'^^^''){F^,Fn 
= eV(E° •B°)(E-B) + --- . 



(18) 



The effective vacuum action arising from the quark box 
diagram in Fig. [2] is thus: 



(QECD) 
eff.O 
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ifzZfl! (E''.B'^)(E-B). (19) 



f f 



It describes the contribution of a heavy quark flavor 
{mf ^ Aqcd ) to the pseudoscalar effective gluon-photon 
coupling. 

The finite temperature effective QED action has been 
studied repeatedly fs^ - ISGj . A detailed account of the 
weak-field expansion of the thermal effective action can 
be found in [35|- The expression for the mixed QED- 
QCD pseudoscalar contribution differs from the pure 
QED result in the same way as described above for the 
vacuum case. The explicit expression, after a Wick rota- 
tion in the integration variable s, is 
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■cff,T 
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Trc,F(^') 



(20) 



The integral is easily carried out, and we obtain as final 
result for the coefficient k defined in (|14p in the high 
temperature limit (T 3> m) : 



«(T) 



651 C(5) V- 
720 (7rT)4 ^ 



217C(5) 
360 (7rT)4' 



(21) 



where the last form includes the contributions from the 
three light quark flavors. We note that the expression 
is independent of the quark masses as long as these are 
small compared to the temperature. 



For the three light quark flavors, the interaction be- 
tween gluons and photons in the pseudoscalar channel 
is intermediated by the iso-singlet pseudoscalar mesons 
?7 and 77' (see Fig. [3]). The effective interactions are 
determined by the electromagnetic and chromodynamic 
anomaly, respectively [STj : 



r 



r 



£. 



r , 



a / cos6' 2 

a f 2 sin 
— — ^ cos 



— ^\/3sin6'r/(E'' •B'^) 

-^V3cos0?/(E'^ -B"^). 
27r/,, 



sin6l) 77(E-B); (22) 

^'(E-B); (23) 
(24) 



V6 



(25) 



Here d w —20° is the flavor singlet-octet mixing angle and 
w 1.2/^ « 157 MeV [111. The effective pseudoscalar 
photon-gluon interaction is then, in the low-energy limit, 
given by: 



£, 



(QECD) 
off.O 



1 



tan d / tan B 

2V2 ml V2\/2 



tan^^ 



^ aa. cos^^^ (E" .B°)(E.B) (26) 



corresponding to the zero-temperature coefficient 
1.46 



(27) 



A useful form that interpolates between the low- and 
high-temperature limits (P71 [^T|) is: 



RiT) 



Ko k{T) 



(28) 



V4 + <Ty 

The temperature dependence of k(T) is shown in Fig. 



C. Hadronic Current 

We first consider the hadronic current in more detail, 
starting from the general form of the "anomalous" inter- 
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FIG. 4: Temperature dependence of the coefficient k in the 
effective pseudoscalar QED-QCD interaction (|14|l . 




FIG. 5: Quark triangle diagram for the electromagnetic cur- 
rent j induced by the presence of a magnetic field B in a 
neutral pion. The cross at the upper corner of the triangle in- 
dicates the insertion of the electromagnetic current operator 

6/7"- 



into a neutral vector meson, e.g. the p'^-meson, which is 
polarized in the direction of the magnetic field and in- 
duces a polarization current among the charged pions in 
the medium. 



action ((T0| . Invoking the general form of the anoma- 
lous current, we obtain for the vector current density in 
vectorial notation: 



Jan 



(dtV) B + (VV) X E] 



(29) 



Identifying V ~ {a/T:f„)(j}, we obtain for the anomalous 
current induced by a magnetic field acting on the pion 
field: 



(30) 



where the dot, as usual, indicates a time derivative. 
What may seem strange about this result is that a neu- 
tral meson field (the neutral pion field) combines with 
the electromagnetic field to generate an electric current. 
However, the effect is easily understood when one recog- 
nizes that an external electromagnetic field deflects the 
quark and antiquark constituents of the pion in oppo- 
site directions causing the neutral pion to become in- 
ternally polarized. The current thus describes the 
polarization current. Before proceeding we note that the 
magnetic field generated by two heavy nuclei colliding at 
relativistic energies is extremely large, and its strength, 
as measured by the product eB, can be of order (see 
Fig. Al in [U). 



The polarization current is highly localized and 
confined to the interior of the neutral pion. It is this cur- 
rent that generates the electric field vector of one of the 
two photons in the decay of the neutral pion. However, 
in order to generate a macroscopic current that can be 
observed as a fluctuating charge asymmetry with respect 
to the reaction plane, we need an effective coupling to 
charged hadrons that can move freely through the nu- 
clear flreball. The relevant Feynman diagram is shown 
in Fig. ini Here the magnetic field excites the neutral pion 
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FIG. 6: Feynman diagram for the electromagnetic current j 
induced by the conversion of a neutral pion into a neutral p- 
meson in the presence of a magnetic field B. In the vector 
dominance model, the electromagnetic current is propor- 
tional to the p-meson field, see Eq. (|32p . 



The effective ir'^V-f coupling, where V stands for any 
neutral vector meson, has been studied extensively |38| . 
Experimentally, the strengths of the couplings are known 
from the partial decay widths of the vector mesons: 
^uj^^o^ = 0.79 MeV, Fpo^^o^ = 0.077 MeV These 
observations imply that there exists an effective interac- 
tion of the form !!4(| 



Zmy 



(31) 



where (f> denotes the neutral pion field and the neutral 
vector meson field {V = u!,p'^)- The effective coupling 
constant gp^^j = 0.58 [i^]. For later purposes, we also 
note the strength of the prj'-f coupling, gp^'^ w 1.31 p^ . 
which can be deduced from the measured decay width 
Try^p^ = ag^y^pUml, « 60 ± 5 keV [4l| . 

Later we will evaluate the fluctuation of the magnet- 
ically induced electric current {ji{x)jk{x')) in a thermal 
hadron gas. The calculation is much simplified by noting 
that the vector meson dominance (VMD) model allows us 
to replace the hadronic electromagnetic current operator 



jf^ by the vector meson operator: 



J' 



9p 



(32) 



For the sake of simplicity, we focus on the p-meson con- 
tribution, because only the p-meson couples to the two- 
pion channel, and pions are by far the most abundant 
hadrons in baryon-symmetric, thermal hadronic matter. 
Using the effective interaction ((3T|) 



C 



cff 



-(B-p)^, 



(33) 



we obtain 

{ji{x)jk{x')) 



Hp 



dy dy' 



P mn 

X {Pt{x)p,n{y)) {Pn{y')Pk{x')) 

X {4>{y)i{y'))B,^{y)B^{y'). 



(34) 



Since we are interested in the matter contribution to 
the in-medium p-meson propagator (p'Hx) p" (y)) , which 
was calculated by Gale and Kapusta [421, '^^ write the 
position-space propagator as 



{p''ix)p-iy)) 



(2^ 



^k{x-y)Di^'^(^k). (35) 



The momentum-space propagator (in the A = oo gauge) 
can be separated into its longitudinal and transverse com- 
ponents: 



fc2 



nL(fc) 



fc2 



nT(fc) 



m?pk'^ 



(36) 



nL/T(fc) denotes the longitudinal and transverse compo- 
nents of the p-mcson self-energy. Since we are primarily 
interested in slowly varying contributions to the induced 
current, i.e. in the limit (k^ , k) ^ nip, the matter contri- 
bution to (p6)) can be approximated as 



-^p,mat(^) 



n 



(mat] 



(fc) 



n 



(mat) 



(fc) 



(37) 
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FIG. 7: Feynman diagram for the electromagnetic current 
j induced by two gluons in the pseudoscalar channel in the 
presence of an external magnetic field B. 



takes the following form for the pseudoscalar two-gluon 
coupling: 



\[aa) 
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Kaas9t(E° •B°)B. 



(38) 



As discussed in Sect. II. A, the pseudoscalar coupling 
of two gluons to the light quark sector in the vacuum is 
dominated by the rj' (or 77) meson. An effective interac- 
tion analogous to (|5T|) then allows the rj' meson to convert 
into a p^ meson in the presence of an external magnetic 
field. The induced anomalous electric current is then 
simply given by the VMD relation ([5^ . The Feynman 
diagram describing this process is shown in Fig. [51 





FIG. 8: Effective hadronic Feynman diagram for the anoma- 
lous electric current j induced by two gluons in the pseu- 
doscalar channel in the presence of a magnetic field B. 



III. ELECTRIC CHARGE ASYMMETRY 



Explicit expressions for the thermal pion gas contribu- 
tions to the p- meson self-energy can be found in [42 1. 



Scenarios for the Reaction Plane Charge 
Asymmetry 



D. Gluon Induced Current 

The diagram in Fig. [71 on the other hand, represents 
a process whereby two gluons induce a current in the 
medium by coupling to a quark- ant iquark pair in the 
presence of the strong magnetic field. Again invoking 
the general form of the anomalous electromagnetic cur- 
rent pip , the anomalous current due to a magnetic field 



We now discuss five different scenarios that may con- 
tribute to the creation of an event-by-evcnt reaction plane 
charge asymmetry in heavy ion collisions due to the in- 
teraction of the magnetic field with the pseudoscalar sec- 
tor of the QCD matter. An overall guidance principle 
for assessing their relative importance is that the cur- 
rent fluctuations are proportional to the square of the 
magnetic field B acting on the strongly interacting par- 
ticles. As analyzed by Kharzeev et al. [2l|, the coherent 
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magnetic field generated by the moving charge density of 
the colliding nuclei in the center-of-mass reference frame 
is strongly peaked at the initial moment of impact and 
then falls off rapidly with an approximate power-law tail 
|B(r)| oc T^^. As a consequence, scenarios occurring 
early in the reaction are strongly favored to be fertile 
environments for the creation of a charge asymmetry. 

Color Glass Condensate ( CCC) Scenario: The mag- 
netic field strength argument suggests that we should 
first and foremost consider the contribution to • B° 
from gluons contained in the saturated gluon wave func- 
tions {i.e. the color-glass condensate) of the colliding 
nuclei. These wave functions contain gluons of all col- 
ors and helicities and thus provide a bountiful supply of 
initial states containin g a nonzero density of topologi- 
cal charge density ntop Indeed, Shuryak and Zahed 
[4^ have argued that the interaction of the colliding sat- 
urated gluon fields will create color field configuration 
of the sphaleron type, i.e. configurations carrying half- 
integer winding number, which then decay into multiple 
quark-antiquark pairs. 

Quark-Cluon Plasma (QGP) Scenario: As already 
pointed out, the operator E" • B° is related to the Chern- 
Simons number density pes of the nonabelian gauge field. 
In thermal equilibrium, pes is given by the thermal wind- 
ing number fluctuations in the quark-gluon plasma. This 
is the process considered by the authors of refs. (2ll[2^. 

Glasma Scenario: Contributions to the anomalous cur- 
rent can also arise from topological charge density fluc- 
tuations occurring during the pre-equilibrium "glasma" 
phase. The magnitude of these fluctuations has been es- 
timated by Kharzeev et al. [i^. When glasma-phasc 
quark pair production [46l |47| occurs in the presence of 
a strong, oriented magnetic field, the produced pairs will 
carry a nonzero electric current in the direction of the 
magnetic field. The generation of an anomalous current 
due to quark-pair production by a chromo-electric flux 
tube was recently investigated by Fukushima et al. [i^ 
in the framework of the effective QCD+QED Lagrangian. 
This process also occurs at very early times, of the order 
of , and thus shortly after the peak in the magnetic 
field strength. 

Corona Scenario: The surface region of the nuclear 
reaction zone never reaches the energy density required 
to form a quark-gluon plasma and thus remains in the 
hadronic phase of QCD matter throughout the reaction. 
Fluctuations of the charge asymmetry in this region will 
be produced by tt — p conversion in the magnetic field. 
This process can act on pious produced during the initial 
impact or even on virtual pious contained in the initial 
nuclear wave functions. 

Hadron Gas Scenario: After re-hadronization of the 
quark-gluon plasma until final freeze-out, the hot QCD 
matter proceeds through a thermal hadronic gas phase, 
where the magnetic tt — p conversion mechanism can op- 
erate throughout the bulk of the fireball. This scenario is 
disfavored by the lateness of its occurrence, as well as by 
its relatively short duration due to the rapid transverse 



expansion of the fireball at late times. 



B. Kinematic Considerations 




A2 Al 



FIG. 9: Schematic view of the collision plane of the two heavy 
ions and the choice of our coordinate system. 



We first outline the strategy for the calculation of the 
charge asymmetry AQ with respect to the collision plane. 
We assume that the collision axis is the x-axis and the 
collision plane the x ~ y-plane, the axis perpendicular 
to the reaction plane is the z-direction. The strength of 
the magnetic field generated by the fast moving nuclei 
acting on the region of space-time causally connected to 
the mid-rapidity region is highly peaked around the col- 
lision moment r = 0. This means that the current j^jyx) 
separating charges with respect to the collision plane is 
concentrated at early times, when the collective trans- 
verse fiow of the matter is small. Wc can thus neglect (in 
first approximation) any collective fiow during the period 
when the charge asymmetry is created. If there were no 
fiow later on, the charge asymmetry would subsequently 
be erased by diffusion. Hence, the collective transverse 
fiow which transports the separated charge further away 
f ! rom the collision plane is an essential component of the 
mechanism. The final charge asymmetry of all hadrons 
emitted above and below the collision plane will be de- 
termined by the charge distribution on the freeze-out hy- 
persurface Ej. We will assume that the local net charge 
density p(x) of the matter at a given point x € E/ can 
be parametrized by an effective local chemical potential 
Pq(x) for electric charge. Finally, it is important to keep 
in mind that the average charge asymmetry (AQ) with 
respect to the collision plane is zero, and only the fluc- 
tuation ((AQ)^) of this quantity can be measured by an 
event average. 

The experiment measures fluctuations of the charge 
asymmetry with respect to the collision plane in momen- 
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turn space, and the chiral magnetic effect creates a fluc- 
tuating cfiarge density asymmetry with, respect to the 
coUision place in position space. In order to convert the 
latter into the former, it is necessary for a correlation 
to exist between the location of the emission point of 
a hadron and its angular distribution. Such a correla- 
tion can have two origins: collective transverse flow or 
the spatial orientation of the freczc-out hypcrsurface. In 
general, one expects both mechanisms to contribute; only 
for isochronous freeze-out or in the absence of collective 
transverse flow would a single mechanism exist. Below, 
we will estimate the amount of expected charge asym- 
metry with respect to the collision plane in these two 
extreme cases, which should give a lower bound on the 
real magnitude of the effect. 

We also note that there are two space-time regions 
of the reaction volume that must be treated separately. 
One is the core region, where the initial energy den- 
sity is high enough to create matter in the deconfined 
phase. Here the chiral magnetic effect induces a cur- 
rent spike, which leads to local charge separation as dis- 
cussed above, followed by transport through the medium 
up to the final freeze-out surface. The other region is the 
"corona" , where the energy density never exceeds the 
critical energy density required to create a quark-gluon 
plasma. Here the evolution proceeds entirely through 
the hadronic phase, and the axial vector current is al- 
ways carried by hadronic interactions. Since there is 
anomaly matching [52[ between the hadronic and quark- 
giuon phases of QCD, the anomalous contribution to the 
induced electromagnetic current remains unaffected, but 
the axial vector correlator will have a different represen- 
tation, in terms of hadronic states. 

Due to the Lorentz contraction, the magnetic field gen- 
erated by the colliding nuclei in the central rapidity re- 
gion is sharply peaked around the collision moment t = 
and falls off by about two orders of magnitude during 
the first fm/c [2lj. This is the time where most of the 
charge separation will be caused by the chiral magnetic 
effect. During this period, the transverse collective flow 
is small, and we will therefore neglect collective flow dur- 
ing the charge separation phase. We then just need to 
calculate the net charge density created within the hot 
QCD medium during the period when the magnetic field 
is strong. We will later discuss how the charge density 
created in the separation process is transported by the 
collective flow field and eventually freezes out into de- 
tected particles. 



C. Asymmetry Generated at Freeze-out 

1. General formulation 

In order to calculate fluctuations in the charge asym- 
metry between the particles emitted in the half-spaces 
below and above the reaction place, we need to know the 
net charge fluctuations on the freeze-out surface. To be 



precise, what we need is the charge asymmetry of the 
particle distribution functions on the freeze-out surface. 
We can then use the Cooper- Frye formula to calculate the 
charge asymmetry of the emitted particles. If we denote 
the freeze-out surface by S/ , the rate of emitted particles 
with momentum p is given by 



The total number of particles emitted is 



N : 



/ (icr^n^(x) 



where 



d^p^Mic,p)0{<j,y 



(39) 



(40) 



(41) 



is the outward directed particle current on the freeze-out 
surface. The electric current and emission of net charge 
is obtained by multiplying each species with its electric 
charge e^. The observable is the difference between the 
total charge emitted into the half-spaces above and below 
the collision plane. In the limit of a completely transpar- 
ent medium, this observable is given by 

AQ = / d^pY.e~sgn{p,) (42) 

i 

= jd^j^ ^^^e./.(x,p)sgn(p,)0(a^p^). 

This quantity is related to the local charge density on the 
freeze-out hypcrsurface S^, which can be expressed as 

P{x)\^f = j d'pY,^^f^i^^^V>)k^■ (43) 



2. Isochronous freeze-out with flow 

To obtain a first result and perform a specific calcu- 
lation, we now specialize to the case of an isochronous 
freeze-out with transverse flow v(x). This means that 
the entire matter volume freezes at a time r = r/, and we 
ignore any opacity for the produced matter. Thus, fluid 
cells in the upper hemisphere, which have a flow vector 
pointing in the upward direction, will be converted to 
an isotropic thermal distribution of particles in the rest 
frame of the fiuid cell. As there is no opacity, many of 
these particles will move in the negative direction per- 
pendicular to the reaction plane. As a result, a large 
part of the effect of charge separation will be washed 
out. Obviously, this represents an underestimate of the 
true magnitude of the effect. 

For isochronous freeze-out we have p^da^ ~ Ed^x and 
9{(Tf^p^) = 9{E) ~ 1. We thus obtain on the freeze-out 
hypcrsurface: 



^Q^ j d^x j d^pY^eifi{x,p;Tf)sgn{p,). 



(44) 
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In order to relate p and AQ, we assume that the particle 
distributions at freeze-out are given by thermal distribu- 
tions with flow and a local chemical potential of electric 
charge, //q(x). For simplicity, we further assume that the 
thermal distributions are well approximated by Boltz- 
mann distributions, and the baryon chemical potential 
vanishes: 

/,(x, p; Tf) = exp [-u,,{^)p^/Tf + e,MQ(x)/T/] , (45) 

where Tf is the freeze-out temperature. Since /iQ(x) is 
assumed to be small, we can expand to first order in this 
quantity. We also find that in the absence of a baryon 
chemical potential, particle and antiparticle distributions 
are equal, and thus the terms independent of /.tg sum to 
zero in both, p and AQ. The expressions linear in pg 
are: 

p(x,T/) =a*q(x) / d3p^e2/(°)(x,p;T/), (46) 
•' i 

and 

AQ = y d'xpQi^) J d3p^e2/(°)(x,p;T/)sgn(p,), 

(47) 

where the superscript "(0)" indicates that we have set 
Pq = 0. We can now eliminate /iQ(x) and obtain an 



expression of AQ in terms of p: 

Ag = / d^xp{x,Tf) (sgn(p^)). 



where 



(sgn(p2)) = 



/'^'pE.ef/f°^(x,p;r/)sgn(p,) 



/d3pE,ef/f^(x,p;r/) 



(48) 



(49) 



We note that the expression on the right-hand-side de- 
pends on position x only through the flow velocity v(x). 
We also note that p is the divergence of a vector, as given 
by eq. below. Thus, the spatial integral over p van- 
ishes, and AQ = in the absence of transverse flow, as 
remarked earlier. 

There are two ways of proceeding from here. The first 
approach assumes that the flow velocity is not large and 
proceeds by expanding the momentum distribution /j'"' 
up to linear order in v: 

/f )(x,p;r/) « exp(-i?p/ry)(l - v(x) -p/T/). (50) 

Here we neglected all terms of higher order in v. Since the 
integrand in the numerator of the square bracket in (|48p 
contains the factor sg'n{pz), only the term proportional 
to Vz contributes. On the other hand, the term v • p docs 
not contribute to the integral in the denominator. We 
thus obtain 



AQw / (fxp{x,Tf) 



!d^pE^ehM-Ep/Tf)v,\p.\/Tf 



m|) 



Jd^pj:.e^cxp{-Ep/Tf) 
Assuming that the sum over particle species is dominated by pions, and using m-n /Tf sa 1, we find 
3 
2 



2E«eMif2(™«/T/)- 



AQ 



d?x p{-K,Tf)vz{yi). 



(51) 



(52) 



The integral does not vanish, because both /o(x) and Vz{^) are antisymmetric with respect to the collision plane. 

In the second approach, which has the advantage that we do not need to make an assumption about the magnitude 
of the fiow, we start by using (j56p to replace p with the anomalous current jan in (j48p . Since both, the numerator 
and denominator of the term in square brackets involve an integral over all momenta, we can make a Lorcntz boost 
into the local matter rest frame at each location. This replaces the particle distributions with those in the rest frame, 
/,;(x, p;Ty ) — > cxp(— E^p/T/), the momentum integrals acquire a Lorentz factor 7(x) with 7 = (1 — v^)^^/^, and the 
argument of the sign function transforms to ^^(x) = Pz + {l — ^)vz'^ • P + ^VzEp, where p'^ denotes the ^-component 
of the momentum in the laboratory frame and p^ is the z-component of the momentum in the comoving frame. We 
thus obtain the following expression for AQ: 

-Jd^pj:^ e\ eM-Ep/Tf) sgn[pi(x)] " 



AQ 



dt / d-\TV •jan(x,t) 



(53) 



jd^pY.,e}cM~EplTf) 

The X dependence of the term in brackets is contained solely in the argument of the sign function. We now integrate 
by parts and obtain 



AQ = 2/ / d3xjan-Vp;(x 



S d^pY.AeM-Ep/Tf)5[p'z{^)] 



Jd^pj:.e^eM-EjTj) 
I 



(54) 



The integral in the numerator is restricted to those par- 



ticles which are at rest in the direction perpendicular to 
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the collision plane. The result ((54)) has a simple inter- 
pretation: It counts the number of charged particles that 
are caused by the anomalous current to change direction 
from downward motion with respect to the collision plane 
to upward motion. The factor 2 accounts for the fact 
that any change in direction causes a gain in the number 
of upward moving charges and a loss in the downward 
moving charges. It is also evident from ([M)) that the gra- 
dient of the flow velocity selects the component of the 
anomalous current that is parallel to it. Thus, if there 
is an anomalous current in any direction other than per- 
pendicular to the reaction plane it will contribute to not 
only the charge separation perpendicular to the reaction 
plane but also that parallel to the reaction plane. 

In the following, we will use the form ([52]) . The fluc- 
tuation of the up-down charge asymmetry is then simply 
given by 

((Ag)2) ^ IJd'xJ dV(p(x,r/)p(x',r/)) 

XV,{x,Tf)v^{x\Tf). (55) 

We now must relate the correlator of the separated charge 
density p to the correlator of the anomalous electric cur- 
rent. 

As discussed above, we can neglect collective trans- 
verse flow during the phase of the collision when most of 
the charge separation with respect to the collision plane 
occurs. The induced charge density p then satisfies the 
continuity equation p + "V ■ jan — 0. We can thus cal- 
culate the charge density by integrating the continuity 
equation: 

p(x,t) = - / diV-ja„(x,<). (56) 
Jo 

The charge fluctuation strength ([55]) can thus be ex- 
pressed as: 

((AQ)2) « ^ y" ' dtdt' J (fx J Sx' (57) 

XVW4X,T/) • (ja„(x,t)ja„(x',t')) ' V'w4x',r/). 

We now introduce the local integrated current fluctuation 
strength tensor C^^: 

C'\^,Tf) = dtdt'J d^x' (j:„(x,i)jt(x',i'))- (58) 

This allows us to express the fluctuations of the charge 
asymmetry as 

((AQ)2)«^ j d'xC'\^,Tf)V,V,{^,Tf)VkV,{^,Tf). 

(59) 

The continuity equation used to obtain the relation 
([55|) does not include the effects of advection and dif- 
fusion. The actual current is composed of the anomaly 
induced current part and the advcctive part; 

j=jan+PV. (60) 



Charge conservation implies 

|^+vV/9=-V-ja„. (61) 

The effect of diffusion on the charge distribution can be 
included by adding the diffusion term (see e.g. in 
the matter frame: 

-DchV^p. (62) 

In any other frame, in which the matter is moving with 
four- velocity u^, the diffusion term has the form 

D,i,{g>''' ~u^u-')d^d,. (63) 

These equations lay the ground for a future comprehen- 
sive and quantitative study of the formation of reaction 
plane charge asymmetry fluctuations by the action of the 
magnetic field on the strongly interacting matter. 

3. Geometric approximation 

The derivation of the fluctuations of charge asymmetry 
leading up to Eq. ([55)) was carried out in the Dch — > 
limit. We will not pursue the consideration of a flnite 
-Dch in much greater detail in this article. We will only 
consider the opposite limit of a very large diffusion co- 
efficient Dch oo, which renders the medium virtually 
opaque. In this limit, the computation of the produced 
charge asymmetry is greatly simplified. 

One may assume that any charge excess in the up- 
per hemisphere of the fireball will contribute to the final 
charge asymmetry with respect to the reaction plane in 
momentum space, and that the same holds for any charge 
excess in the lower hemisphere, i.e., the large diffusion co- 
efficient does not allow the charges to move back to the 
opposite hemisphere. This assumption clearly overesti- 
mates the effect, because not all charges in the upper 
hemisphere will eventually be emitted with an upward 
component of the momentum. However, it is useful to 
obtain an upper limit on the size of the expected effect 
which, together with the result of the calculation assum- 
ing an isochronous freeze-out, can help bracket the pre- 
diction of a more sophisticated calculation. 

In this limit we define the charge asymmetry fluctua- 
tion geometrically as 

{{^Qf)^-( <fx ( dV(p(x)p(x')). (64) 

J z<0 J z>0 

Using the continuity equation, as before, we can rewrite 
this as a double integral over the charge current through 
the reaction plane: 

((Ag)2)= j d^xd^x'5{x^)5{x'^){f(x)f{x')). (65) 

Note that the minus has disappeared, because one of the 
currents flows upward through the reaction plane; the 
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other one downward. As already noted, in conjunction 
with the Dch calculation this geometric estimate 
allows us to bracket the uncertainty in our results arising 
from the flow and diffusion dynamics of the final state. 
Of course, this docs not provide for an estimate of the 
inherent uncertainty in the approximations leading up to 
the calculation of the local anomalous current density. 




B 



FIG. 10: Feynman diagram for the electromagnetic current- 
current correlator (j(2;)j(a;')) generated by gluon fusion in the 
pseudoscalar channel in the presence of an external magnetic 
field B. 



IV. ESTIMATE OF THE CHARGE 
ASYMMETRY FLUCTUATIONS 



In this section, we derive semi-quantitative estimates 
for the charge asymmetry generated by the chiral mag- 
netic effect in collisions of two Au nuclei at the top RHIC 
energy. We first consider the contribution from gluon fu- 
sion as depicted in Fig. [71 where the two gluons are part 
of the incoming nuclear gluon distributions. The second 
contribution corresponds to the case where the gluons are 
part of the thermal distribution in the deconfined phase. 
Finally, we consider contributions from the hadronic pro- 
cess depicted in Fig. [B] in a thermal environment with a 
temperature T sa 160 MeV. 

We begin with the calculation of the completely trans- 
parent medium, i.e. the Dch — ^ limit, followed by the 
geometric estimate for completely opaque (Dch ^ oo) 
medium for both partonic processes. We only derive the 
estimate for the upper bound of an opaque medium for 
the hadronic contribution. We finally note that one is 
really interested in the quantity 



((A(iV+-iV_)f) 



((AQ)^ 



(66) 



because experimentally one measures the difference in the 
number of equally or oppositely charged particles emitted 
into the two hemispheres, not the charge difference. 



A. CGC Scenario 

1. Anomalous current correlator 

We first consider the contribution depicted in Fig. |S1 
which describes the anomalous electric current generated 
by the fusion of two gluons from the colliding nuclei, 
in the pseudoscalar channel, due to the presence of the 
strong magnetic field carried by the colliding nuclei. We 
now estimate the integrated magnitude of this current 
and compare it with the overall quark pair multiplicity. 



We start from the following expression for the current- 
current correlation function, shown graphically in Fig. 1101 



^9 prj' ■ 



cosf 



9p 



X J d^yd^wd^w'd^y' {p^{x)pHy)) B,{y) 

X {ii{y)i{w)){{^'^ ■B'^){w){^' ■^'){w')) 

X {r^\w')r^'{y'))B,Xy'){p^{y')p\x')). (67) 

This expression describes the process, whereby a gluon 
pair fuses into a virtual ry'-meson, which is converted into 
a neutral p-mcson by the intermediation of an external 
magnetic field. Equation (|67p assumes the vector meson 
dominance relation ([32]) for the electromagnetic current. 

One might argue that it would be more appropriate 
to describe this process at the partonic level, because 
the final state into which the current is imbedded is a 
pre-equilibrium quark-gluon plasma and not a hadronic 
environment. However, a first estimate is more easily ob- 
tained using the effective hadronic representation of this 
process discussed in Sect. II. A. We expect that the mag- 
nitude of the hadronic and partonic mechanisms shown 
in Figs. [5] and [3] is related by the parton-hadron duality 
property of QCD [s^ . A confirmation of this conjecture 
by an explicit evaluation of the analogous partonic dia- 
gram would be desirable; we will not attempt it here. 

We now note that all elementary correlators in (|67p . 
being controlled by hadronic states with masses of the 
order of 1 GeV, are highly localized compared with the 
size of the nuclear interaction region. It thus makes sense 
to approximate them by space-time delta functions, in 
particular: 



{p:,{x)p^{y)) 



^S^'Hx-y), 
—5^^\y~w). 



(68) 
(69) 



We also note that the magnetic field Bz{y) is localized 
in the longitudinal direction due to Lorentz contraction 
of the nuclear charge distribution to a region of charac- 
teristic width i?/7 « 0.05 fm (for a Au nucleus at top 
RHIC energy of 100 GeV/u). Since this is even less than 
the characteristic width of the gluon cloud of the color 
glass condensate (approximately « 0.1 fm), we can 
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average the magnetic field over a region of width in 
the longitudinal direction. Using the result derived for 
Bint derived in the Appendix, this yields 



ei?2 « QseBj 



int 



2Za 



67 



(70) 



where b denotes the impact parameter of the nuclear col- 
lision. We note that the average value Bz is nearly inde- 
pendent of the transverse coordinates within the nuclear 
reaction volume. 

Introducing the abbreviation 



3p^/^ 12{Zayai cos^ 6> 6^72 



52 (2^/,)2m2m2 i?6 



(71) 



we then obtain the following expression for the current- 
current correlator: 

0.(x)j.(x')) « C{{E^ ■ B^){x){E' ■ B'){x')}. (72) 

We now have to evaluate the correlator of the 
color glass condensate fields in the two colliding nu- 
clei. Because the color fields in a fast nucleus are of 
the Weizsacker- Williams kind, i.e. approximately plane 
waves with orthogonal polarizations of the chromoelectric 
and chromomagnetic field strengths, the relevant contri- 
butions to the pseudoscalar invariant E° • B" arise when 
E° and B° originate in different nuclei: 

= ^([E^(:r).B^(x)][E^(x')-B^(x')]) 

= E(^""(^)4o(^')> {Bl0{x)Blp{x')), (73) 

where the indices a, /3 = 1, 2 count the two colliding nu- 
clei and i,j denote the spatial vector indices. 

We can relate these matrix elements to the gluon dis- 
tribution function G'(^) in the colliding nuclei as follows. 
Color and rotational symmetry imply that for a single, 
fast moving nucleus the matrix element is diagonal in the 
indices (a, 6) and («,j), where i,j are directions trans- 
verse to the beam, and that the chromo-electric and - 
magnetic correlators are equal: 

(Kr(x)£;^(x')> = {Bnx)B'^^ix')) 

_ SijSab 



2(7V2 _ 1) 

^ijSab 
4(7V2 - 1) 

Si j Sab 

m? 1) 



{Et{x)Et(x')) 
m{x)Et{x')+Bt{x)Bnx')) 



{F''+\x)F''+\x')), 



(74) 



where = {Ef + e^^jB^) / \/2 is the transverse chro- 

modynamic field strength in the direction along the light 



cone. The gluon distribution function in a proton is de- 
fined as a light-cone Fourier transform of precisely this 
matrix element: 



eG(e) 



e-«?'^"" {p\F-+\u-)F-+^mp). (75) 



Neglecting nuclear modifications of the gluon distribution 
in the proton, the nuclear matrix element (|74p is related 
to the matrix element in the proton as: 



{F''+'{x)F''+''{x'))a 



{p\F^+^(u)F^+^iO)\p)P^^ 



, (76) 



where x = {x + x')/2 and u ~ x — x', and /o(x±, x^) is the 
nuclear density distribution, which is normalized to A. 
Nuclear modifications can be taken into account by using 
a modified gluon distribution function of the nucleon. 
Inverting the Fourier transform in (|75p we obtain the 
relation: 

{F^+'{x)F''+%x'))a 
= p+ j d^e-'^iP^^^'-^"HG{0p{5c^,x-). (77) 

Since we are here interested in the charge asymmetry at 
mid-rapidity for heavy-ion collisions in the RHIC energy 
domain, the relevant value of the Bjorken parameter ^ of 
the gluon distribution function is ^0 ~' 10"^ at a scale of 
order ^ m^, « 1 GeV^. Approximating 

6G(Co) « 3 (78) 

as a constant in this kinematic range, the Fourier integral 
evaluates to 



p+ j d^e-'^P^^'' ) = 2T:Six- - x'- 



(79) 



If we also approximate the Lorentz contracted nuclear 
density distribution as 



p(x_L,.T ) K Ta{'x.\_) 6{x ), 



(80) 



we obtain the desired connection between the nuclear ma- 
trix element of the gauge field strength fluctuation: 

{F°'+\x)F''+'{x'))a 

= 2iT5{x-)5{x'~)[£,oG{io)]TA{5^,_). (81) 

Combining this result with the analogous expression 
for the second nucleus traveling in the negative light-cone 
direction, we finally obtain for the matrix element (|73l) : 

{[E^{x)-B%xm\x')-B\x')]) 
(27r)2 



m - 1 



5{x-)5{x'-)5{x+)5{x'+) 

X KoG(eo)]'T^A(x±;b), (82) 
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where 



TAA(x±;b) 



Ta X 



^ ) . (83) 



A similar expression is obtained when one evaluates the 
gauge field strength correlator (173]) in the color glass con- 
densate model. Following Lappi l49j, one obtains 



^KoGlCo)] 



(84) 



where Qs(Co) is the nuclear saturation scale and R is the 
nuclear radius. The values ([75]) for [^oG(^o)], = 0.3, 
and i? = 7 fm for a Au nucleus correspond to a choice of 
the saturation scale, = 1.7 GeV^. 



2. Isochronous freeze-out approximation 

We now make use of the isochronous freeze-out result 
([57]) in the I?cii limit for the final charge asymmetry 
fluctuation. In doing so, we apply the approximation [B • 
Vvz] ~ [BzdzVz], which is exact for the type of transverse 
flow field we will be considering. The current-current 
correlator from Eq. ([7^ can then be used to yield the 
charge fluctuation: 
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e^C / d'^xd'^x'[dzVz{x)][d'zVz{x')] 



((E" •B'')(x)(E'' •B'')(.t')). 



((AQ)^ 



(85) 

We remind the reader that y, z denote the transverse 
directions with respect to the beam in the four-vector 
x^^ — {t,x,y,z). Making use of the expression ([5^ for 
the color field correlator, we can perform the integrals 
over t,t',x,x' and write ([55]) as ((AQ)^) = {9/4)e^CI 
with 



x[eoG(^o)]'rAA(x±;b) 



J][9>.(xl)] 



(86) 

The double integral can be expressed as one over xj^ and 
one over the difference Ax^ = x^ — x^. The integral 
over Ax_L requires some explanation. At first sight, the 
integrand depends on Ax^ only very weakly through the 
collective flow field (x^). It turns out that this is an ar- 
tifact of our use of the integrated gluon distribution func- 
tion G(^), where we expressed the nuclear matrix element 
([77]) in terms of the gluon distribution. For x^ ^ x'j_ 
we should have used the unintcgrated gluon distribution 
G(^, k^) instead. This would have led to the replacement 



[eoG(Co)] 



KG(e,ki)] 



?kj_ - Ax 



(87) 



Assuming that the transverse fall-off of the nuclear un- 
intcgrated gluon distribution is controlled by the satu- 
ration scale, this substitution could be approximated by 
the analytic expression 



KoG(eo)] KoG(eo)]e" 



(88) 



Neglecting the variation of the flow velocity gradient over 
distances of order 1/Qs, the integral over Ax^ can now 
be carried out, and we obtain: 



[N^ - 1)(2Q,)2 



(89) 



To proceed further, we need to make explicit assump- 
tions for the transverse flow profllc and for the nuclear 
density proflle. For the transverse flow velocity field at 
freeze-out we assume a self-similar, linear profllc of the 
form 



v_l(x^) = w/xx/i?. 



(90) 



where R is the nuclear radius. This choice implies a con- 
stant gradient: dzVz = Vf /R. Again for simplicity, we as- 
sume the nuclei to be approximately homogeneous, solid 
spheres with density p and radius R, which gives us 



Ta{^±) ^ 2pJ R^ - X 



(91) 



for the nuclear thickness function. A rough approxima- 
tion for the nuclear overlap integral is found as: 



dyTAA{v;b) 



with 



f{h) 
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b 



fib) 



Putting everything together, we obtain 



lib) 



97rA^[^oG{Co)]hj 
8(iV2 - l)(Q,i?2)- 



f{b). 



(92) 



(93) 



(94) 



Making use of the correspondence ((84|). the analogous 
expression in the color glass condensate model is 



CGC) 



(b) 



9{Nl 



l)v) 



(8^)3a2 



Qlfib)- 



(95) 



We note that (jM]) contains a factor 1/a^, which ac- 
counts for the nonperturbative gluon density of the color 
glass condensate. This factor cancels against the factor 
Ofg in the result for the matrix element (|7ip . The final 
result for the charged particle number asymmetry fiuc- 
tuation, based on expression (l94l) is thus 



((A(7V, 



N-))') = -^CI{b) 

_ 243 TT gl^,^ {Za)^alcos^e 
= ~8~ gl (27r/,)2m2,m2 



(iV2 _ 1)Q2^8 ^2 



IT, fib). 



(96) 
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Alternatively, using the color glass condensate model ex- 
pression (|95p . the result reads: 



243 g'prj'-y {Zafcos^e 



(87r)3 g2 (2^/^)2^2,^2 

Inserting the numerical values gpr^'-y = 1.31, Qp = 5.03, 
= 1.3 GeV, 7 = 100, Z = 79, i? = 7 fm, we obtain: 

((A(7V+ - N^)f) « 5 X 10-^ v) ^ f{h). (98) 

3. Geometric approximation 

For comparison we now calculate the charge asymme- 
try in the geometric model in the D^h — > oo limit. Ac- 
cording to our reasoning at the beginning of this section, 
leading to the expression (|65p for the charge asymmetry 
fluctuation, we need to evaluate the expression 



j d^xd^x' 5{z)5{z'){{¥.'' ■-B''){x){¥}' ■B,''){x')), (99) 



where we remind the reader that y, z denote the trans- 
verse directions with respect to the beam in the four- 
vector = {t,x,y,z). Making use of the expression 
([5^ for the color field correlator, we can write (jHH) in 
the form 



(2^)^ 
m-1 



dydy'[^oG{^o)?TAA{mb). (100) 



The integral over Ay = y — y' requires some explana- 
tion. At first sight, the integrand does not depend on 
Ay. This is an artifact of our use of the integrated gluon 
distribution function G(^) where we expressed the nu- 
clear matrix element (|77p in terms of the gluon distribu- 
tion. For y ^ y' we should have used the unintegrated 
gluon distribution G(^, k^^), which would have led to the 
replacement 



KoG(^o)] 



dkydk, [eG(e,ki)] e 



(101) 



Assuming that the transverse fall-off of the nuclear un- 
integrated gluon distribution is controlled by the satu- 
ration scale, this could be approximated by the analytic 
expression 



KoG(Co)]^KoG(^o)]e-^='^-^'l. 



(102) 



The integral over Ay can now be carried out, and we 
obtain: 



(27r)^KoG(eo)]^ 



dyTAA{y\b). 



(103) 



Modeling the nuclei as homogeneous, solid spheres with 
density p and radius _R, we have 



TA{^i_) = 2pjm -^i 



(104) 



for the nuclear thickness function and a good approxima- 
tion for the nuclear overlap integral is: 



dyTAA{y;b) 



with 



Putting everything together, we obtain 
6A2KoG(eo)]' 



m = 



(105) 



(106) 



(107) 



Making use of the correspondence ([Mj) . the analogous 
expression in the color glass condensate model is 



(108) 



We note that (|108p contains a factor l/ofj, which ac- 
counts for the nonperturbative gluon density of the color 
glass condensate. This factor cancels against the factor 
a1 in the result for the matrix element (|7ip . The final 
result for the charge asymmetry is thus 

{{A{N+-N.)f) = cm 

^ ^^g%'^{Zaf al cos'' e 



9p (27r/,,)2m2,m2 

QsR'iN^ - 1) 

Alternatively, using the color glass condensate model ex- 
pression (jlOSp . the result is independent of the strong 
coupling constant: 

{{A{N+-N^)f) = CI^^^^Hb) 
9 glrj'-y {Zaf cos' e 

g^p (271/^)2^2,^2 

x{N^-lh'^^f{b). (110) 

Inserting the numerical values gpri'-y = 1.31, gp = 5.03, 
Qs = 1.3 GeV, Z = 79, R = 7 fm, we obtain: 



{{AiN+ N^)f) « 1.7 X 10-3 f{b). (Ill) 



The result (|llip obtained in the geometric model is 
larger than the result (|98p found for the isochronous 
freeze-out by a factor ~ 35/^2. This difference arises 
from an additional factor vi/RQs in the freeze-out model 
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and constitutes a generic effect of using the flow profile 
function versus the geometric separation, which should 
occur independently of the specific mechanism underly- 
ing the anomalous current. The same factor separating 
the lower and upper bound of the predicted effect should 
thus be present, as well, for the other charge separation 
mechanisms considered below. 

We note that the expressions (P7|) and (|110p depend 
sensitively on the masses of the intermediate hadronic 
states, rrijf and nip. If these are strongly modified by the 
collision environment on a time-scale of order 0.2 fm/c, 
which is likely the case, the magnitude of our estimate 
will be affected correspondingly. In order to treat such 
modifications realistically, a more microscopic approach 
will be needed. 



B. QGP Scenario 

We next consider the quark-gluon plasma scenario, in 
which fluctuations of the winding number density of the 
gauge field are driven by thermally assisted transitions 
across the so-called sphaleron barrier [5^ between vacua 
characterized by neighboring integer winding numbers. 
The thermal diffusion rate per unit volume for the Chern- 
Simons number in QCD, also called the "strong sphaleron 
rate" , was numerically computed by Moore [59| with the 
result: 



1 



- sph 



VAt 

„2 



d xd x' {pcs{x)pcs{x')) 



^ J d^x ((E» • B°)(.t)(E'' • B^)(0)) 

(112) 



a 
(2^) 

100 a! T'' 



with an estimated uncertainty of a factor of two. In the 
expression above V denotes the volume and At the in- 
tegration time. In order to extract the winding number 
density correlator from this relation, we assume that the 
spatial and temporal correlation lengths are much shorter 
than the size of the volume occupied by the quark-gluon 
plasma and its lifetime. This assumption is well justi- 
fied, as the characteristic size of winding number carrying 
fluctuations at high temperature is given by the magnetic 
length scale of order (g^T)~^, and the characteristic time 
of a saddle point transition is of order {g'^T)~^ , as dis- 
cussed by Arnold et al. [g^I- All of these scales are of 
order 1 fm or less. We thus express (|112[) in the form 

{{E--B''){x){E''-B''){x'))^S\x-x') — r,p^. (113) 

Starting from Eq. psp we again make use of our previ- 
ous result (|57p for the final charge asymmetry fluctuation: 



((AQ)= 



9 



d X d x' (naas) 



(114) 



X [B{x)-\7v,{x)][B{x')-Vv,{x')] 
X dtdt'{{-E'' ■B-){x){-E'' -B'Xx')). 



We now partially integrate with respect to t and t' and 
insert our explicit expression (jll3p for the pseudoscalar 
gluon density correlator: 



((AQ)^ 



(37ra)^ 



d'^x 



sph 



[Bix) -Vv.ix)] 



(115) 



As we found in Section II, k oc in the high- 

temperature phase of QCD. which implies that the prod- 
uct (K^Fsph) increases like T"^ as the temperature falls. 

We now make the assumption that a reasonable esti- 
mate of the product [B • Vvz] is obtained by retaining 
only the z-component of the magnetic field: [BzdzVz]- 
The time behavior of at late times at the center of 
mass of the colliding nuclei was investigated in ref. [2l| 
(see Fig. A2). For proper times t > tq — 0.5 fm/c the 
results can be approximately parametrized as 



eB{T) ^ eBo 



s{b) 



(116) 



where eBo w 300 MeV^ and the exponent s depends 
somewhat on the impact parameter 6, ranging from s ~ 2 
for 6 = 4 fm to s « 3 for 5 = 12 fm. We further assume 
that the temperature changes with proper time as in the 
boost-invariant longitudinal expansion model: 



(117) 



with initial temperature Tq = 400 MeV. Finally, we as- 
sume that To and eBg are roughly constant over the 
transverse area; this leads clearly to an overestimate, 
which can easily be improved if so desired. We now in- 
sert these expressions into (jll5l) and introduce coming 
space-time variables (t, r/, xj^): 



drj 



^0 Jr„ ^ T 



2s+i 



(118) 



where t/j denotes the hadronization time of the quark- 
gluon plasma. We note that v^^ according to ([57]) . is to 
be evaluated at the freeze-out time Tf. Dividing both 
sides by and inserting numbers we obtain 



dmN+-N^)f 

dj] 



2 X 10" 



2,s-i 



X / d^xi^ (dzVzf 



(119) 



where we have taken the upper limit of the integral to 
infinity. 

To proceed further, we again assume that the trans- 
verse fiow velocity field at freeze-out can be approximated 
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by the hnear self-similar linear profile ((90|) . This choice 
implies dzV^ = Vf/R and thus 



(120) 



Further neglecting the constant (—1/3) in the denomi- 
nator, we obtain out final result for the charged particle 
number asymmetry fluctuations with respect to the re- 
action plane in the QGP scenario: 



dmN+-N^)f) 

drj 



3 X 10 



(121) 



We finally note that when one evaluates the CGC 
scenario in the geometric approximation, one obtains 
an ill-defined product of three delta functions in the z- 
direction: 6{z)5{z')6{z — z'). This means that one can- 
not neglect the correlation length of the Chern- Simons 
number density in the z-direction, but must assign it a 
nonzero value A ^ (p^T)"^. The delta function product 
then becomes 5{z)5{z') / \, and the final result acquires a 



factor (-R/A) instead of the factor 



C. Corona Scenario 

We finally consider the corona scenario, where the 
anomalous current is generated by the hadronic process 
shown in diagram Although we will, for definite- 

ness, evaluate the current-current correlator in a thermal 
ensemble, we do not insist that this is an excellent ap- 
proximation at the relevant early times when the mag- 
netic field achieves is maximal strength. Nevertheless, we 
hope that the result obtained in this way gives an upper 
estimate of the magnitude of the expected effect. Denot- 
ing the location of the ptttt vertex by y and that of the 
pBir^ vertex by z, the current operator is given by 



fix) = i.gp,,($*(a;)a^$(x)-(a^$*(a;))$(x)) 

{<i>*{y)d''<i>{y) - (9''<i>*(2/))$(2/))p.(2/) 



X / d 



X ^e.pxs J d''zp\z)d'cj,{z)FP%z). (122) 

We recall that <() denote the charged and neutral pion 
fields, respectively. 

Making us of the translation invariance for the thermal 
current-current correlator, we can simplify the expression 
(I65p for the charge asymmetry fluctuation to read 



((AQ)2) = j d^xd: 



x'6{x^)S{x''){f{x)fm, (123) 



which will serve as the starting point of our calculation. 
Translation invariance also makes it convenient to con- 
sider the Fourier transform j'^{k) of the current. In the 
thermal ensemble, the Fourier transform of the current- 
current correlator is given by the "<" unordered corre- 
lator, because we are interested in the effect due to an 



incoming neutral pion and not in vacuum fluctuations. 
We define: 



n^^(fc)= / d''xe-''^f{x)fm 



(124) 



We note that only the longitudinal component of the cur- 
rent vector, and thus only the longitudinal component of 
the p- meson field contribute to the matrix element, be- 
cause we have used the continuity equation to relate the 
charge density to the divergence of the current. Further- 
more, only the component F^^ = of the magnetic 
field perpendicular to the reaction plane contributes to 
the correlator. For moment, we neglect the time- and 
space-dependence of the magnetic field B. We also note 
that the integration over the coordinates (x and y) span- 
ning the reaction plane reduce the momentum range of 
interest to = ky = and = E^^^ (k) = y/k"^ + . 

After some algebra, one then obtains the following ex- 
pression for n'^'^(fc): 



n^-^(fc) 



2 ^ n2 '^''T'Trik'^) 



n(fc)2 



{k^-ml-n,{k)Y 

where Ilp{k) is the p-mcson self-energy and 



(125) 



m) = 



d^p 2(2p3 - k^){p°k^ - p^k°) 



(27r)4 (po)2-(i?(±)(p))2 

n(jP) 

[{po-EAp)r-{Ei^\p-k))^] 



(126) 



Here E'i^'' (p) denotes the (in-medium) on-shell energy of 
a charged pion, while i?,r(k) denotes the on-shell energy 
of the neutral pion initiating the current. Note that the 
integral over p^ in the last equation accounts for the sum 
over Matsubara frequencies. The integral over the four 
poles of the integrand in (|126p yields two pairs of identical 
residues, leaving the expression: 

2p3 _ fc3 



{E{k)-E{p)y -E{p-k)^ 
2p^ + k^ 
{E{k)+E{p)f -E{p + k)\ 



(127) 



where we have dropped the sub- and superscripts on the 
pion on-shell energies to reduce the cluttering of the equa- 
tion. Since the only non-zero component of k is k^ , we 
can go to polar coordinates in the integral over p and 
integrate over the azimuthal angle. This leaves a two- 
dimensional integral over p = \p\ and the polar angle 
cos 6. This integral simplifies considerably if we neglect 
the pion mass relative to the temperature T. Since the 
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main contribution to the anomalous current comes at 
early times, this may provide for a reasonable approx- 
imation. The final result is: 



n(fc) 



1 

1 



dpp 



6 



dcos 9 — niv) 

-1 P 



(128) 



where we have written = for clarity. 

We now need to substitute (|128p into Eq. (|125p and 
transform back to coordinate space. However, doing so, 
we need to take into account the time dependence of the 
magnetic field. The full expression to be evaluated is: 



((AQ)'> - 



9p-K7T 



7l2 (m2 
dt dx dy 



dt' 



2^ 



ikt 



X K{k)B^{x+'^)B^{x^'^). (129) 

We will use the following approximation for the magnetic 
field generated by a single nucleus (see (|A8[) in the ultra- 
relativistic limit v ^ \): 



Zebj e{R~y)e{R~ z) 



47r (^2 +^2(^^i)2) 



3/2 



, (130) 



where the — (+) sign is for a right-moving (left-moving) 
nucleus. The magnetic field of both nuclei is given by 

We begin with the Fourier integral, which is approxi- 
mately of the form 



Ji(fc) 



^ikt 



dt'- 



TT e 



-kR/j 



kR 

7 



kR 

7 



(131) 



Next we perform the integration over the current momen- 
tum k. Noting that -//R > T for 7 = 100 and T ss 150 
MeV, implying that kR/j <C 1 where the remainder of 
the integrand is large, the integral simplifies to: 



dk k 
27r 2 
3 



1 



64i?5 



dk' 







h{k) 

k _ n^T'- 
- 1 ~ 1287^' 



(132) 



Next, we perform the integrals over t,x,y. Since the 
magnetic field depends on 7(x =F t), the nuclear volume 
is Lorentz contracted in at least one light-cone direction; 
a reasonable estimate may thus be 



dt dx dy 



i?3 



(133) 



Collecting all factors, we get: 

((AQ)-). '-"~'' (^)°^. (134, 

Inserting Z = 79, gp^^ = 0.6, Qp^^ = 5.92, T sa 150 MeV, 
we finally obtain 



mN+-N^)f) = 



3x10-5—, (135) 



which is approximately 50 times smaller than our geomet- 
ric estimate (|llip for the contribution due to the color 
glass condensate scenario. We note once more that our 
result likely constitutes an overestimate of the true mag- 
nitude of the corona contribution, because a thermal pion 
gas will not be formed until some time (maybe 1 fm/c) 
after the onset of the collision, when the magnetic field 
strength has already subsided significantly. We did not 
take this delay time into consideration in our estimate. 



V. SUMMARY 

In this work, we have explored possible sources of elec- 
tric charge asymmetry fluctuations caused by the inter- 
action of the coUisional magnetic field with the highly 
excited QCD matter created in a relativistic heavy ion 
collision. We first analyzed the various microscopic mech- 
anisms that contribute to the anomalous current and 
showed that the chiral magnetic effect, i.e. charge sep- 
aration along the direction of the magnetic effect, oc- 
curs in both, the partonic and the hadronic, phases of 
QCD. In the partonic phase after thermalization, the 
anomalous current is dominantly produced by winding 
number carrying gauge field configurations interacting 
with the magnetic field via a thermal quark loop. Be- 
fore thermalization, the winding number carrying gauge 
field configurations can interact with the magnetic field 
via a virtual eta' (or rf) meson, which electromagneti- 
cally converts into a p-meson. In the thermal hadronic 
phase, the anomalous current is predominantly gener- 
ated by the electromagnetic n — p conversion process 
7r° + 7 ^ p'' ^ 7r+ + TT-. 

It is evident from our analysis that, contrary to pre- 
vious claims, reaction plane charge asymmetry fluctu- 
ations do not require "local" or global parity violation 
for their formation, as the constituent processes of these 
various mechanisms - except the thermal winding num- 
ber fluctuations in the quark-gluon plasma phase - are 
well known from ordinary hadronic physics. We discussed 
five distinct scenarios for the creation of a reaction plane 
charge asymmetry fiuctuation, viz. the color glass con- 
densate (CGC) scenario, the quark-gluon plasma (QGP) 
scenario, the glasma scenario, the corona scenario, and 
the hadron gas (HQ) scenario. We analyzed two of these, 
the CGC and the corona scenario, in some detail and 
found that the estimated magnitude of the expected ef- 
fect in both scenarios is much smaller than the effect 



18 



observed by the STAR experiment. This suggests that 
the observations may be due to some other process. 

Rigorous predictions of the reaction plane charge 
asymmetry fluctuations will need to include in-medium 
effects on the interactions involved in the various mecha- 
nisms, such as modifications of the meson masses and ef- 
fective couplings. They will also need to track the charge 
density transport from the creation by the anomalous 
current through the expanding matter up to the freeze- 
out hypersurface. We showed how the final charge asym- 
metry fluctuation can be calculated from the fluctuation 
of the charge density at freeze-out. 

While it is difficult to quantitatively ascertain the rel- 
ative magnitude of the influence of these effects, it is rea- 
sonable to expect that in-medium modifications of meson 
propagators may somewhat enhance the charge asym- 
metry, whereas dissipative transport mechanisms active 
during the expansion phase will likely suppress the asym- 
metry. It is instructive to attempt a very rough estimate 
of upper bound by which in-medium modifications of the 
masses of the intermediate states involved in the charge 
asymmetry creation could enhance the expected effect. 
As an example, we consider the CGC scenario. Here the 
magnitude of the charge asymmetry fluctuation is con- 
trolled by the coefficient C, eq. (|7ip . which involves the 
factor (27r/^m^' nip)'^ in the denominator. These factors 
represent scales relating to the properties of intermediate 
quark-anti-quark states. In the normal QCD vacuum, 
these will be limited below by the QCD confinements 
scale Aqcd ~ 200 MeV; in a thermal medium, they will 
be limited below by the thermal mass scale of order gT, 
which is also of the order of 200 MeV or more under con- 
ditions reached at RHIC. An optimistic upper limit to 
the possible enhancement due to in-medium modification 
of the effective mass scale of quark-antiquark excitations 
would thus be (27r/^m^/ mp)V(200 MeV)^ « 10"^. 

It is thus not inconceivable that QCD processes occur- 
ring in a high energy density environment could enhance 
the charge asymmetry fluctuations by several orders of 
magnitude above our estimates. In order to result in 
such an extreme enhancement of the charge asymme- 
try fluctuations, however, the processes responsible for 
it must take hold on a time scale less than 1 fm/c so that 
they can act during the period of peak or near-peak mag- 
netic fields. Furthermore, if the chiral magnetic effect is 
experimentally confirmed as the source of the observed 
event-by-event fluctuations, it thus would imply the ex- 
perimental observation of a form of highly excited QCD 
matter in which the creation of local winding number 
fluctuations is strongly enhanced compared to the nor- 
mal QCD vacuum. 
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Appendix A: Magnetic field estimate 

The magnetic field of the two colliding nuclei is calcu- 
lated from the retarded vector potential 



A{x) 



D^ctix - x')jNix') 

D,.ct{x - x')pn[x' - vt')-fv, 



(Al) 



where v is the velocity of the colliding nucleus, 7 the 
Lorentz factor, and pN the static nuclear charge density. 
Here we have neglected any "slow-down" of the valence 
quarks as the two nuclei collide. In a more detailed calcu- 
lation one would want to use a more detailed treatment of 
the dynamics of the matter in the fragmentation regions 
of the colliding nuclei. 

Suppose a point charge e is moving along the x-axis 
and its x-coordinate is given by a; = vt. See Fig. 12] for 
a schematic view of the kinematics and the coordinate 
system. The magnetic field B at r = (x, y, z) is given 
by 



B = X r' 



(A2) 



where 



(7(3- - vt), y, z). 



(A3) 



When a point charge e is moving along y — -^^ 
magnetic field created by the charge is given by 



±f, the 
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tV X r 



±1 



where 



{"f{x-vt), j/Tt:, z). 



(A4) 



(A5) 



In the following, we set a; = and consider the magnetic 
field as a function of y and z. The z-component of the 
magnetic field is the superposition of the ^-components 
of the magnetic fields created by the two point charges: 

B,^B^+^+Bi-\ (A6) 



where 



ev"f 



3/2 ■ 



(A7) 



For nuclei, we adopt the approximation that the charge 
Ze is distributed uniformly with a radius R. Then the 
above formula (|A7|) is replaced by 



Zavj (2 T 2/) 



{r± < R), 



Zavj{lTy)^ {r±>R), 



(A8) 
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where 
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y T 2 ) + 2^ + "^7^t^ 



1/2 



(A9) 
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Next we calculate the time integral of 



dteB^{x ^0,y,z,t) (AlO) 



and note that the result is independent of the product jv. 
As typical values applying to a mid-eentral Au+Au col- 
lision at RHIC, we take 5 = 7 fm, R ~ 7 fm, and Z = 79. 
The time-integrated field eBi^t is shown in Fig. [TT] as 
function of the position z perpendicular to the reaction 
plane for y ~ (top panel), and as function of the po- 
sition y within the reaction plane for z ~ 0. As the 
figure shows, the time-integrated magnetic field is nearly 
independent of the location within the overlap region of 
the two nuclei. The impact parameter dependence of the 
time-integrated magnetic field at the center {y = z = 0) 
of the transverse plane is shown in Fig. [121 in comparison 
with the approximate formula 



eBu 



2.32Za— r. 



(All) 



FIG. 11: Time-integrated magnetic field strength eBint{y, z — 
0) (top panel) and eBint{y ~ 0, z) (bottom panel) in units of 
MeV for the parameters b — 7 fm, R — 7 fm, and Z = 79. 
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